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Abstract 

We derive the exact solution of the Boltzmann kinetic equation for 
the three-dimensional Lorentz model in the presence of a constant and 
uniform magnetic field. The velocity distribution of the electrons reduces 
exponentially fast to its spherically symmetric component. In the long 
time hydrodynamic limit there remains only the diffusion process governed 
by an anisotropic diffusion tensor. The systematic way of building the 
Chapman- Enskog solutions is described. 
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1 Introduction 

The model proposed by Lorentz in 1905 to describe the electrical conduc- 
tivity in metals has become an inexhaustible source of results concerning the 
foundations of the kinetic theory. Within the model the electrons are considered 
as small hard spheres with electric charge (— e) propagating among randomly 
distributed immobile hard spheres representing the atoms. The electrons do not 
interact with each other. The elastic electron-atom scattering is at the origin of 
the electric resistivity of the medium. 

An excellent review of the subject up to 1974 can be found in 2 . Since 
then the research keeps bringing new results. More recent applications of the 
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model, in particular in the theory of dynamical systems have been described in 
[3]. One of the relatively recent discoveries was the non-Mar kovian character of 
the evolution in two-dimensions in the presence of a magnetic field perpendic- 
ular to the plane of motion. Even in the Grad limit where the fraction of the 
volume occupied by the scattering atoms is very low the kinetic equation keeps 
containing memory effects and does not reduce to the Boltzmann form 

We consider here the three-dimensionsional dynamics where the mean free 
path of the electron A is inversely proportional to the number density of the 
scatterers n and to the scattering cross section na^ with a equal to the sum of 
the atomic and electronic radii. Keeping A fixed while the volume fraction na^ 
approaches zero defines the Grad limit 

A = l/nira^ ~ const , , 

In the absence of external fields the Boltzmann equation applies in the regime 
Moreover, it can be rigorously solved. This fact has been exploited in to 
compare the evolution predicted by the Boltzmann equation for general initial 
conditions with that of the special class of Chapman-Enskog solutions involving 
exclusively hydrodynamic modes. 

Our object is to generalize the analysis of |S] to the case where the electrons 
are acted upon by a constant and uniform magnetic field B. As we consider 
the motion in three-dimensions the non-Markovian effects studied in |4j| do not 
occur and the Boltzmann equation yields an adequate description in the low 
density limit (^) . The electrons move between collisions under the action of the 
Lorentz force F giving rise to the acceleration 

F 

— = a; X V = cjB x v (2) 

m 

where B — i?B, |B| = 1, m is the electronic mass, and ui — eB/m is the 
cyclotron frequency. 

We denote by /(r, v, i) the probability density for finding the electron at 
time t at point r moving with velocity v. The Boltzmann equation applied to 
the Lorentz model (also called Boltzmann-Lorentz equation) reads 

fir, V, i) = - { [P/](r, V, t) - /(r, v, t) } (3) 
r 

The collision frequency \/t ~ v/\i& constant under the assumed dynamics. In- 
deed, the speed v = |v| of the electron does not change during the free cyclotron 
motion and is also conserved by the elastic collisions. The gain term in the right 
hand side of ((SJ represents the effect of isotropic scattering which averages the 
distribution / over all directions leaving the spherically symmetric part 

[P/](r,t,,t) = j ^/(r,v,t) (4) 

V = v/v represents a point on a unit sphere in the velocity space so that c?v is 
the solid angle measure. 

It turns out that the initial value problem for the kinetic equation ^ can be 
explicitly solved. We begin by deriving in section 2 the evolution of the velocity 
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distribution which permits to calculate the diffusion tensor via a Green-Kubo 
formula. The complete exact solution for the distribution / is derived in section 
3. Section 4 is devoted to the discussion of the hydrodynamic diffusion process 
governing the evolution in the long time limit. In section 5 we indicate the 
method of a systematic construction of Chapman-Enskog solutions whose time 
dependence is entirely defined by the diffusive mode at any time scale. The final 
section 6 contains discussion and comments. 



2 Velocity distribution and diffusion tensor 

2.1 Approach to spherical symmetry 

The evolution equation for the velocity distribution 

(/.(v,t) = J dr/(r,v,t) 
follows by integrating Q over the position space. One finds 



d , ,9 



(v,t) = --[0(v,t)- [P<^](t;,t)- 

T 



Applying the projection 



to both sides of © yields the relation 

d_ 
dt 



(5) 



(6) 



(7) 



The spherical part [Pc/i] of the velocity distribution remains unchanged in the 
course of time. This permits to rewrite (jS)) in the form 



d . . d 

— + X V • — 
ot ov 



X(v,t) =0 



(8) 



with 



X(v, t) = exp(t/T) [ <^(v, t) ~ [Fcp] {v,0)] (9) 

The solution of ^ is straightforward as it consists in propagating x(v, t) back- 
ward in time up to i — along the free cyclotron trajectory 

v(-i) = n^{-ujt) • V = V|| + n^{-ojt) ■ v_L (10) 

where 72.g (a) is the rotation of angle a around the axis B. The velocity compo- 
nent parallel to the magnetic field V|| = (v ■ B)B is not affected by the cyclotron 
rotation whereas the perpendicular component = v — V|| rotates with con- 
stant angular velocity lo. We thus have in general 

X(v,<) -x(v(-<),0) 

which in view of @ yields the solution 

with 

0o(v) =0(v,t = O) 

The velocity distribution becomes spherically symmetric exponentially fast with 
the relaxation time r = X/v. Its projection [P0o](w) does not change in the 
course of time and represents the asymptotic equilibrium state. 



(11) 

(12) 
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2.2 Diffusion tensor from a Green-Kubo formula 

The complete solution of the initial value problem permits to evaluate 
the conditional probability P(v,t|w,0) for finding the electron with velocity 
v(i) = V at time t provided it had velocity v(0) = w at time t = 0. One finds 



S{v — w) 



(13) 



The knowledge of P(v,t|w,0) suffices to determine the diffusion tensor 
through the Green-Kubo formula 



-t-oo 



D^., = / dt (v,,(i)v,(0))e, (14) 
Jo 

where the velocity auto-correlation function is given by 

(Vp(t)v,(0))e9 = / '^^ / '^WV^W,P(v,t|w,O)0eg(?«) (15) 

Here 

0eg(w) = S{w - wo)/47ruo 

denotes the spherically symmetric equilibrium state. 

The (5-distributions greatly facilitate the evaluation of the integral in H15|l 
which is readily reduced to 

(16) 

Here [TZ-^{ujt)^ are the elements of the rotation matrix corresponding to rota- 
tion of angle uit around the B-axis. The calculation of the diffusion tensor (|14|l is 
conveniently performed in the coordinate system whose z-axis is oriented along 
the magnetic field. One then finds the diffusion matrix with diagonal elements 

D± = D^^ = Dyy = ^"o^YT^^ ' ^" " " ^^'^^ 
and the off-diagonal elements 

D.,y = -Dy,^-^VIT,^^^, D,,^Dy,=0 (l8) 

Equations H17 |l . p8|l show the reducing effect of the magnetic field on the diffu- 
sion process perpendicular to the magnetic field. The fact that the correspond- 
ing diffusion coefficient D±^ decays at large values of the field as (o^r)"^ has 
been well established in plasma physics [Hj. The formulae H17|l . IjlSI) have been 
also derived with the use of the appropriate Langevin equation [J] . The physical 
content of expressions H17(l and H18|l is discussed in detail therein. 

Our main object in this paper is to exhibit the time scale separation in the 
evolution toward equilibrium. In contradistinction to the rapid relaxation Ijllfl 
of velocities, there also appears a by far slower hydrodynamic spreading of the 
electrons in the position space governed by the diffusion tensor. 
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3 Solution of the Boltzmann equation 

On the left hand side of the Boltzmann equation Q we find the operator 

-+v.- + (a;xv).-J (19) 

which generates collisionless trajectories in the magnetic field. The backward 
in time motion starting from the phase space point (r,v) reads 

v{-t) = r-V||^-i7^g(|) ■[vi(-t)-vi] 

w{-t) = n^{-ojt) ■ V = v|| + n^{-LJt) ■ (20) 

(see equation (llOH . The operator 1)19(1 vanishes when applied to the trajectory 
(I2U|I which permits to rewrite J^J) in the integral form 

/(r,v,t) = e-*/Vo(r(-t),v(-t)) + - f dt'e- ^[^f]{v{-{t - t')),v,t') 

T Jo 

(21) 

with /o(r, v) = /(r,v,t = 0). Applying to H21() the Fourier transformation 

/(k,v,t)= / dre-* V(r,v,t) (22) 



we find 

/(k,v,i) =cxp{-i/T-ik- [v||i + x(vi,i)]}/o(k,7lB(-cji)-v) (23) 
1 f* 

+ - / dt'exp{-(t-t')/T-ik. [v||(t-t')+x(vi,(t-t'))]}[P/](k,i;,t') 



(24) 



where the shorthand notation 



(v±,^) = ^7^B (I) •[vj.(-0-v±] 



has been introduced. 

The way to the solution consists now in using the fact that the Laplace 
transformation 



/(k,v,z) 



dte--*/(k,v,0 



(25) 



factorizes the convolution in time present in 123(1 . The integral form 1(21(1 of the 
Boltzmann equation for the Fourier-Laplace transform / takes the form 



/(k,v,z) = [L/oJ (k,v,z)+VXk,v,z)[P/](k,«,z) 

where the direct effect of the initial state is contained in 

Jq (k,v, z) = 



(26) 



(27) 



f + 00 



dt exp <j — ( z + i + ik ■ V[| ] t — ik ■ , i) ^ /o (k, T^g (— wi) • v) 
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whereas the function 
■0(k,v, z) = ■ 



dt exp 



z -\ h ik ■ V|| j t — ik • x(v^ , t) 



(28) 



multiphes the spherical part of the transformed distribution /. Owing to this 
factorized structure the apphcation of the projection P to both sides of l|2t)l) 
yields a closed equation for [P/](k, v, z) providing the formula 



[P/](k,v,z) 



PL/o 


(k,'y,z) 


1 - [P'v 


;](k,i;,z) 



(29) 



By inserting the solution H29|) into (|26|l we arrive at the complete solution for 
the Fourier-Laplace transform of the distribution function 



/>,v,z) 



(k,v, z) 



^(k,v, z) 
1 - [Pi/;](k,w,z) 



PL/o 



(k,w,z) 



(30) 



The exact evolution in the course of time can be inferred from H3(J|I by taking 
the inverse Laplace transformation 



dz 

/(k,v,t)= / — exp(zt)/(k,v, z) 
27ri 



(31) 



The contour of integration in H31() is a line parallel to the imaginary axis lying 
to the right of all singularities of the analytic function /(k, v, z). In order to 
analyze the modes of the time evolution we have thus to examine the singularities 
of hmctions L/q, PL/q and ijj together with zeros of the hmction (1 — P?/;) where 



dt exp < — z 



[PV'](k,w,z) = (32) 
— ^ < j- y — exp {— ik ■ V||t — ik • x(v^, f)} 



ith 



x(vi, t) = - (f - ■ v± - (I) ■ 
(see definition (|^ ). 



(33) 



4 Separation of time scales: appearance of the 
hydrodynamic mode of diffusion 

The formulae H27|l . H28|l show that L/o(k, v, z) and "0(1^, v, z) are Laplace trans- 
forms of bounded (periodic) functions evaluated at the point {z + \/t + i\s.- V||) 
whose real part equals (Rez -\-1/t). It follows that all their singularities in the 
complex z-plane are located in the half-plane Re(z) < — 1/t (the same holds for 
PL/o ). The modes of evolution related to these singularities correspond thus 
to a rapid decay dominated by the exponential decrease ~ exp(— t/r). 

A qualitatively different behavior comes from the singularity occuring when 
the denominator in (|29|l vanishes. Consider the equation 

[Pi/'](k,w,z) = 1 (34) 



6 



Its explicit form is conveniently written in the coordinate frame whose z-axis 
is parallel to B and where vector k lies in the plane xOz with coordinates 
(fcj_,0,fc||). We find 

exp - z + - U / / (35) 



T I V T/jJo 2 7o 27r 

1 . ^Lut 



X exp ■j^— i<fc||W cos 9 — i — sin J 2k±v sin 9 sin (f)j = l 

where {v sin 9 cos 0, u sin 9 sin (/>, t; cos 9) are the coordinates of the velocity vector. 
Denoting by ZhydO^) the implicit function defined by 1)35(1 one checks that at 
k = there is a unique solution z^ydO^ = 0) = 0. Upon expanding the integrand 
in 1(35(1 in powers of k± and fc|| we find the asymptotic expansion of z;iyd(k) for 
k — > in the form 

Zhydik) = - [i^||fc2 + D^kl\ + o(fc2) (36) 

where the diffusion coefficients _D|| and D± coincide with those derived from the 
Green-Kubo relation ((14(1 (see ((17(l ). o(fc^) denotes a term vanishing faster than 

The isolated zero a.t z = ZhydO^) of the analytic function {1—Pip) corresponds 
to the hydrodynamic pole in the solution 1(30(1 . Indeed, in the vicinity of 2 = 
and for k ^ one finds the asymptotic representation 

' - ' (37) 



l-FiP{k,v,z) Tiz - Zhydi'k)) 

The evaluation of the inverse Laplace transformation 131(1 will thus yield the 
contribution 

Ay<i(k, v,t) ^ exp[-i(i?||fcjj + D^kl)] (38) 

which represents the slow diffusion process spreading the electrons in the po- 
sition space. Their number density in the Laplace-Fourier representation (at a 
given modulus of velocity) follows directly from ((30(1 by applying the projector 
P 

-('-'-'-^- i-tp'Kk,.,.) ^''^ 

Here the hydrodynamic zero z — z^yd determines entirely the evolution for long 
times i 3> 1/t. 

The magnetic field does not change the kinetic energy of the electrons. 1(39(1 
is thus the density of a conserved quantity which implies the appearance of the 
hydrodynamic mode. The situation here is analogous to that studied in "7 in the 
absence of B where the existence of Chapman-Enskog solutions holding all along 
the evolution has been proved by construction. In the next section we show how 
to construct systematically Chapman-Enskog solutions in the presence of the 
field B. 



5 Chapman-Enskog solutions 

The only hydrodynamic field in the dynamics of the Lorentz model is the number 
density 1(39() . In order to investigate the possibility of the existence of Chapman- 
Enskog solutions depending on time only via n(r,v,t) we follow here the ideas 
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developed in j2| where such solutions have been explicitly constructed. In the 
Fourier representation we thus look for distributions of the form 



/(k,v,t) = [P/](k,i;,t) + [(l-P)/](k,v,i) 
= ■^HKv,t) + L{k,v\h{.,t)) 



(40) 



If / solves the linear Boltzmann equation 

/(k,v,t) = 



— + ik • V + (a; X v) • — 
at av 



/(k,v,i)-[P/](k,^;,i) (41) 



then a/ is another solution for any number a. It follows that L(k, v|n(., t)) is 
a linear functional of h 



ai(k, v|n(., t)) = L(k, v|an(., t)) 
Following Hauge we assume the simplest local linear relation 

L(k,v|n(.,t)) = i^(k,v)n(k,?;,t) 
and look for the Chapman-Enskog solutions of the form 
/(k,v,<) - [l + ^^(k,v)][P/](k,^;,i) 



(42) 
(43) 
(44) 



In order to derive a closed equation for the function F we insert (|44|l into l|41|l . 
We get 



(9 3V 
(1 + F)— + ik • v(l + F) + (a; X v) • — 



In view of the relations 

P[F] = 0, P[ik • v] = and 



[P/](k,t,,t) = --F[P/](k,^;,t) 

(45) 

(46) 



(uj X v) 



dF 
dv 



= 



the application of the projection P to (|45|l yields 



d 



— + [V(ik-v)F]}[Pf](k,v,t)^0 



which determines the evolution of the spherical part of the distribution 
[rf]{k,v,t) = [P/o](k,«)exp{-i[P(zk.v)i^](k,w)} 



(47) 



(48) 



The time evolution of the complete distribution H44|l has been thus also deter- 
mined 



/(k,v,t) = [l + F(k,v)][P/o](k,t;)exp{-t[P(ik.v)i^](k,i;)} 



(49) 



When inserted into l|45|l the equation H47|l permits to derive the nonlinear equa- 
tion determining the function F 



1 , , d 
- + {uj xv) ■ — 
T av 



F = {[P(ik • ^r)F] - ik • v} (1 + F) 



(50) 
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Notice that by putting k = in (|47|l we find 

^ + lim[P(ik-v)F]i [P(/)](t>,i) = (51) 

ut k— >0 J 

where 0(v, t) is the velocity distribution. H51|l should coincide with the previ- 
ously derived condition lO expressing the invariance of the spherically symmet- 
ric component of 0. We conclude that 

lim[P(ik- v)i^l = (52) 

Taking then the limit k in H49() we find that the velocity distribution 
does not change in the course of time which is possible only for spherically 
symmetric initial conditions (see This property of the Chapman-Enskog 

solution imposes on F the condition 

F(k = O,v)=0 (53) 

Indeed, as the relation (|44l) must hold at any time, in particular at t = 

/(k,v,0) = [l + F(k,v)][P/](k,i;,0), (54) 

and the restriction of the class of initial conditions to those whose velocity 
distributions are spherically symmetric enforces (|53(l . 

Coming back to the Laplace-Fourier transform / we get from (|49f) the for- 
mula 

/(k, V, z) = [1 + ^^(k, v)] [P/o] (k, v) ^ , (55) 

z + [F(iK ■ •v)}' \ (k, V) 

The Chapman-Enskog distribution /(k, v, z) has thus a single hydrodynamic 
pole at 

z = ZhydCk) = - [Fiik ■ v)F] (k, v) (56) 

which determines entirely the evolution in the course of time. 

A subtle point should be noted here. Owing to the presence of spatial gradi- 
ents the Chapman-Enskog distribution is not spherically symmetric as function 
of V all along the evolution towards equilibrium. However, as a consequence of 
()53|l its k = component possesses the spherical symmetry at any moment of 
time. 

The power series representation of the function F can be constructed in a 
systematic way in terms of convenient variables (see (ESI) 

si = \i± ■ -v^ = k±vsm6cos(j) (57) 

52 = (k^ X v^) • B — k±v sm6 sin(j) 

53 — k|| • V|| — k\\vcos 9 

We write 

F = F^-^^ + F^"^^ + . . . (58) 

where is the n-th order polynomial in variables (si, 32,33). The action 

of the operator C = {uj x v) ■ d/dv appearing in (|50|l on these variables is 
particularly simple 

Csi = —UJS2, Cs2 = wsi and £33 = (59) 
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which greatly facihtates the calculations. The nonlinear equation H5U|I implies 
the chain of linear equations 



1 , . d 
- + a; X V ■ — 



1 , . d 

- + (a; X v) • — 
T av 



F(i)(k,t)) = -ik- V (60) 
(k, -y) = [¥{ik ■ v)i^(i)] (k, -y) - ik • vF^^^ (61) 



which can be consecutively solved. For the first two terms of the expansion one 
finds 

F(i)(k,v) = -ir S3 + — — — (si+cjTSa) 

1 + W^T^ 



(62) 



F(2)(k,v) 



2 „2 



r s 



(63) 



(1+Cj2r2)2 



[2siS3 + tjr(3 + Li;^T^)s2S3] 



(l + a;2^2)(i+4^^2^2 



[(1 + tj2T^)Si + 3tJTSlS2 + 3cjV2s2] 



According to H56|l the second order term in the expansion of the hydrody- 
naniic pole is given by 



(2) 
^hyd 



»(ik- v)F(^'](k,w) 



(64) 



which reproduces the already derived formula H36|l . For symmetry reasons any 
term f (2«) H5;j|) is a polynomial of even order in the variables si, S2 and S3. 
Therefore, by virtue of definitions (|57|l . i^^^n) jg invariant by the change v — v 
and [P(ik • v)F'^2n)](k^ -y) — : the small k expansion of the hydrodynamic pole 
contains only terms of even orders in k. 

Notice that is in fact the equation for the number density 



d_ 

dt 



[V{ik ■ v)F] 



n(k, t) = 



(65) 



In the small k limit we can use here the first order approximation F^^K In the 
position space this gives the diffusion equation 



(66) 



A II and Aj_ are the Laplace operators corresponding to space variables conjugate 
to the wave vectors k|| and kx, respectively. 

In principle, one could continue the evaluation of higher order terms going 
beyond the classical diffusion equation H66|l . In fact, the determination of the n- 
th order polynomial F*-"^ in the series H58|l is straightforward once one knows all 
polynomials of lower rank. To prove the existence of Chapman-Enskog solutions 
there remains the open question of convergence of expansion (|5^ 
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6 Discussion and comments 



The exact solution (|30|l of the Boltzmann equation for the three dimensional 
Lorentz model in the presence of a constant and uniform magnetic field allows 
to describe precisely the dynamical evolution of the system much like in the 
previously studied case |S] of free motion between collisions. The picture which 
emerges from our analysis bears strong resemblance to the case with B = 0. 
This is clearly related to the fact that by taking the limit uj = eB/m 
in various predictions following from (|30|l we recover the results of We 
could thus reveal the time scale separation between the rapid exponential decay 
of the deviation from the spherical symmetry in the velocity distribution 
and the slow diffusive spreading in the position space due to the hydrodynamic 
pole H3t)|l . This was achieved by showing that the Laplace- Fourier transform 
/(k, V, z) contained a single pole approaching z = for k — > whereas the 
rest of the singularities of / remained within the half-plane Rez < — l/r. We 
have also checked that the Green-Kubo formula (|14|l gave the exact value for 
the diffusion coefficients D±, by identifying them independently in the small 
wave vector expansion of the hydrodynamic pole 1)36(1 . Of course, there is here a 
new effect compared to The magnetic field introduces an anisotropy into the 
diffusion process absent in the field free case. For the diffusion in the direction 
perpendicular to B the Lorentz model predicts the decay D± ~ (wr)"^ in the 
limit of strong fields ojt 1. The diffusion parallel to B looks the same as in 

There is one point where the presence of the magnetic field makes the prob- 
lem much more difficult than in the case of B = 0. ft concerns the possibility 
of solving the nonlinear equation (|50|l , necessary for the explicit construction of 
the Chapman-Enskog solutions. In the limit oj ^ Q H5()|l takes the form 

-F^{[¥{ik-v)F]-ik-v}{l + F) (67) 

T 

and can be readily solved yielding 

"^("'-)'--°- ...cot(J)+..k.v -^ 

It is quite remarkable that the solution of the nonlinear equation H67|l is so 
easily obtained. It permitted in [5] to provide the proof of the existence of the 
Chapman-Enskog solutions by explicit construction. This was however not the 
case with our equation H5()|l . All we were able to do was to find the way to a 
systematic expansion of i^(k, v) in a power series of the variables k±, fcy. The 
first two terms (|62|l and (|63|) of the expansion have been exhibited in the text. 
They become (—irk • v) and [— (rk • v)^ -I- (rkv)'^ /3], respectively, in the w — > 
limit reproducing the first two terms in the expansion of ((68|l . The still open 
question is the convergence of the complete series ((58(1 . 

Let us finally remark, that the further generalization consisting in including 
also an electric field would require a modification in the dynamics of the Lorentz 
model, as keeping elastic collisions with fixed scatterers excludes the existence 
of a stationary state because of the unbounded absorption of energy . 
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